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Abstract 

We compare the mass spectra and string tensions of SU (2), SU (3) and SU (4) gauge the¬ 
ories in 2-1-1 dimensions. We find that the ratios of masses are, to a first approximation, 
independent of Nc and that the remaining dependence can be accurately reproduced by 
a simple 0(1/A^) correction. This provides us with a prediction of these mass ratios for 
all SU{Nc) theories in 2-1-1 dimensions and demonstrates that these theories are ‘close’ 
to Ac = oo for Ac > 2. We also find that, when expressed in units of the dynamical 
length scale of the theory, the dimensionful coupling is proportional to 1 /Ac at large 
Ac. We confirm that these theories are indeed confining in the limit Ac —> oo. We 
describe preliminary calculations in 3-1-1 dimensions which indicate that the same will 
be true there. 
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1 Introduction. 


The proposal to consider SU(Nc) gauge theories as perturbations in powers of 1/iVc around 
Nc = oo is an old one Q. If one assumes confinement for all then the phenomenology 
of the 517(oo) quark-gluon theory is strikingly similar to that of (the non-baryonic sector 
of) QCD 1^, 1^. This makes it conceivable that the physically interesting SU{3) theory 
could be largely understood by solving the much simpler SU{oo) theory The fact 
that the lattice SU{oo) theory can be re-expressed as a single plaquette theory Q, has 
provided the basis of a number of interesting computational explorations (for a review 
see 1^). Unfortunately this latter scheme makes no statement about the size of even the 
leading corrections to the Nc = oo limit, and so gives us no clue as to how close AT; = 3 is 
to Nc = oo. 

In this paper we calculate the properties of SU{Nc) gauge theories for several values of 
Nc and explicitly determine how the physics varies as Nc increases. We note that there 
are good reasons for believing that the inclusion of quarks would not alter any of our 
conclusions (except in some obvious ways). We consider both 2+1 and 3-1-1 dimensions. 
The former calculations are much the more precise and it is there that we will be able 
to make some firm statements. In T) = 3 -|- 1 our conclusions will be similar but more 
tentative. 

While one might naively expect that the D = 2 + 1 and D = 3 + 1 gauge theories 
would be so different as to make a unified treatment misleading, this is not in fact so. 
Theoretically the D = 2 + 1 theory shares with its 17 = 3 -T 1 homologues three central 
properties. Firstly, at short distances, the dimensionless coupling becomes weak in both 
cases. In the D = 2 + 1 theory the coupling, has dimensions of mass so that the effective 
dimensionless expansion parameter on a scale a will be ag^ which vanishes linearly with 
distance (the theory is super-renormalisable). In the D = 3 + 1 theory the coupling 
vanishes logarithmically with distance (asymptotic freedom). Secondly, at large distances 
both theories appear to be confining, with a non-perturbative linear potential between 
fundamental sources. Thirdly it is the value of the coupling that sets the overall mass 
scale in both cases. In 77 = 2-|-l this arises directly because has dimensions of mass. In 
77 = 3 -|- 1 it does so through the phenomenon of dimensional transmutation: the classical 
scale invariance is anomalous, the coupling runs and this introduces a mass scale through 
the rate at which it runs (i.e. the parameter). In addition to these general theoretical 
similarities, the calculated spectra also show some striking similarities. All this motivates 
us to believe that a unified treatment makes sense. 

The 77 = 2 -|-1 analysis is based on our calculations over the last few years of the properties 


2 


of SU{Nc) gauge theories with Nc = 2,3 and 4. In D = 3 + 1 what we have done is to 
perform some SU{4:) calculations to supplement what is known about SU{2) and SU{3). 
Our strategy is the very simple one of directly calculating the mass spectra of these 
theories and seeing whether they are approximately independent of Nc- The calculations 
are performed through the Monte Carlo simulation of the corresponding lattice theories. 
In the D = 2 + 1 case the calculations are very accurate and we are able to extrapolate our 
mass ratios to the continuum limit prior to the comparison. In the 3+1 dimensional case 
our SU{4:) calculations are not good enough for that, and our comparisons with SU{2) 
and SU (3) are correspondingly less precise. 

This study was originally motivated by the observation that the C = + sector of the light 
mass spectrum turned out to be quite similar in the D = 2 + 1 SU (2) and SU (3) 0 
theories. (This also appears to be the case in Z) = 3 + 1, although there the comparison 
is weakened by the much larger errors.) If the reason for this is that both are close to the 
Nc = oo limit, then this provides an economical understanding of the spectra of SU{Nc) 
gauge theories for all Nc, i.e. there is a common spectrum with small corrections. 


A second reason for studying Nc ^ oo is that models and theoretical approaches are 
usually simpler in that limit. For example, the flux tube model of glueballs |^, would 
naively appear to be identical for Nc > 2. However, because the model does not incorporate 
the effects of glueball decay, it should in fact be tested against the Nc ^ oo spectrum since 
it is only in that limit that there are no decays. A second example is provided by the 
recent progress in calculating the large Nc mass spectrum using light-front quantisation 


techniques |1C]. 

The calculations in this paper have all been performed with the standard plaquette action 
using standard Monte Carlo techniques. The lattice spacing, a, is varied by changing the 
dimensionless bare inverse coupling, (5, which appears in the lattice action. In the pure 
gauge theory the states are necessarily composed entirely of gluons and we shall therefore 
refer to them as ‘glueballs’. Some of the SU{2) results have been published |^, ]T| and a 
paper on the other work is in preparation. A brief summary of the work in this paper has 


been presented elsewhere \12^\. 


2 2+1 dimensions. 

We begin with the string tension, a, since it turns out to be the most accurate physical 
quantity in our calculations. We use smeared Polyakov loops |Q, to obtain a^a for 
several values of the lattice spacing a. We then extrapolate the lattice results, using the 
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asymptotic relation (3 = ^jag^, to obtain the continuum string tension in units of 



( 1 ) 


The results for SU{2) Q, SU{3) and 5C/(4) Q, in ZD = 2 + 1 are shown in Table ||. We 
immediately see that there is an approximate linear rise with and we find that we can 
obtain a good fit with 



( 2 ) 


We obtain a similar behaviour with the light glueball masses (see below). 

Some observations. 

• For large Nc, eqn(2) tells us that oc g‘^Nc. That is to say, the overall mass scale of 


the theory, call it /r, is proportional to g'^Nc- In other words, in units of the mass scale of 
the theory 



( 3 ) 


While this coincides with the usual expectation based on an analysis of Feynman diagrams, 
we note that here the argument is fully non-perturbative. 

• The string tension is non-zero for all and, in particular, for Nf. ^ oo (when ex¬ 
pressed in units of g'^Nc or the lightest glueball masses - see below). This confirms the 
basic assumption that needs to be made in 4 dimensions in order to extract the usual 
phenomenology of the large-iVc theory ii- 

• In the pure gauge sector one expects (again from an analysis of Feynman diagrams) 
IQ] that the first correction to the large-W limit will be 0(l/iV^) relative to the leading 
term. The fit in eqn(2) is indeed of this form. We note that if we try a fit with a 
0(1/Nc) correction instead (which would be appropriate if we had quarks) then we obtain 
an unacceptably poor (corresponding to a confidence level of only ~ 2% in contrast to 
the ~ 45% we obtain for the quadratic correction). We may regard this as providing some 
non-perturbative support for the diagram-based expectation. 

• The coefficient of the correction term is comparable to that of the leading term, sug¬ 
gesting an expansion in powers of 1/W that is rapidly convergent. Indeed one has to go 
to Nc = 1 before the correction term becomes comparable to the leading term. While the 
SU{1) theory is completely trivial, we note that the 17(1) theory has a zero string tension 
(in the sense that jg^ = 0 in the continuum limit). 

In addition to the string tension we have calculated part of the mass spectrum. In partic¬ 
ular we have calculated the masses of the lightest particles with quantum numbers 
J = 0,1,2, P = ± and C = ±. In some cases we have calculated some of the excited 
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states for given . Note that since we are in Z) = 2 + 1, states of opposite parity are 
degenerate as long as J 7 ^ 0. This degeneracy is broken by lattice spacing and finite vol¬ 
ume corrections. We will present our results separately for the P = + and P = — states 
so as to provide an explicit check on the presence of any such unwanted corrections. 


We begin with the C = + spectrum since the SU ( 2 ) spectrum does not contain C = — 
states. In this case we have masses for three values of W, and so can check how good is a 
fit of the kind in eqn(2). In Fig. [l|. we plot the ratio mc/g'^Nc against 1/N^ for a selection 
of the lightest states, G. On this plot a fit of the form in eqn( 2 ) will be a straight line 
and we show the best such fits. As we can see, the data is consistent with such a 1/N^ 
correction being dominant for W > 2. However what is really striking is the lack of any 
apparent W dependence for the lightest and 2 ++ states. 


In Table ^ we present the results of fitting the C = + states to the form 

_ p . ^slope 


(4) 


where Roo = j^ ■ (Note that the errors on the slope and intercept are highly 

correlated.) For each state we show the confidence level of the fit. These are acceptable 
suggesting once again that for > 2 a moderately sized correction of the form \INI 
is all that is needed. Note that since the variation with Nc is small, the exact form of the 
correction used will not have a large impact on the extrapolation to W = 00 (except in 
estimating the errors). 


These calculations confirm our earlier claim that the physical mass scale at large Nc is 
g'^Nc- So if we consider ratios of niQ to (as was explicitly done in |]l^) we will find that 
they have finite non-zero limits as Nc ^ 00 : that is to say, the large-Ac theory possesses 
linear confinement. 


For the C = — states we only have masses for 2 values of Ac and we cannot therefore 
check whether a fit of the form in eqn(4) is statistically favoured or not. However given 
that such a fit has proved accurate for the C = + masses and for the string tension down 
to Ac = 2 it seems entirely resonable to assume that it will be appropriate for Ac > 3 for 
the C = — masses. Assuming this we obtain the results shown in Table ^ for the Ac = 00 
limit and for the coefficient of the first correction. The ‘lever arm’ on this extrapolation 
is, of course, shorter than for the C = + states and that leads to correspondingly larger 
errors. 


The results in the Tables provide us not only with values for the various mass ratios in 
the limit Ac —> 00 but also, when inserted into eqn( 2 ), predictions for all values of Ac. 

Finally we remark that we have also calculated the deconfining temperature, Tc, for SU{2) 
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| 14| and for SU{3) IQ]. Extrapolating as in eqn(4), we find 


= 0.1745(52)+ 




m 


( 5 ) 


Of course, extrapolating from Nc = 2,3 is less reliable than extrapolating from Nc = 3,4. 


3 3+1 dimensions. 

Our knowledge of 4 dimensional gauge theories is much less precise. As far as continuum 
properties are concerned, quantities that are known with reasonable accuracy include the 
string tension, the lightest scalar and tensor glueballs, the deconfining temperature and 
the topological susceptibility. As in 3 dimensions, the SU (2) and SU (3) values are within 
~ 20% of each other, which encourages us to investigate the 5t/(4) theory so as to see 
whether we are indeed ‘close’ to N^, = oo. Of course S'f7(4) calculations are much slower 
in H = 3 + 1 and the results we present here are of a preliminary nature. 

We use the standard plaquette action, and so our first potential hurdle is the presence 
of the well-known bulk transition that occurs as we increase the inverse bare coupling, 
(3 = 2Nclg^, from strong towards weak coupling. To locate this transition we performed 
a scan on a 10^ lattice and found that it occurred at /3 = 10.4 ±0.1. This corresponds to 
a rather large value of the lattice spacing, a, and so does not lie in the range of couplings 
within which we shall be working, i.e. /3 =10.7,10.9 and 11.1. 

Our calculation consists of 4000,6000 and 3000 sweeps on 10^,12^ and 16^ lattices at 
/3 =10.7,10.9 and 11.1 respectively. Every fifth sweep we calculated correlations of (smeared) 
gluonic loops and from these we extracted the string tension and the masses of the lightest 
0'''+ and 2'^~^ particles, using standard techniques [^]. These are presented in Table Q. 
We also calculated the topological susceptibility, a^Xt =< where L‘^ is the 

number of lattice sites and Q is the total topological charge. (Note that in 11 = 2 ± 1 there 
is no such charge.) The charge Q was obtained using a standard cooling method, just as 
in SU{3) [^. The calculations were performed every 50 sweeps. Overall this corresponds 
to rather small statistics and the errors are therefore unlikely to be very reliable. 

We see from Table Q that the most accurate physical quantity in our calculations is the 
string tension, a. Can we learn from it how varies with N^, just as we did in H = 2± 1? 
We focus on a simple aspect of this question: if we compare different SU {Nc) theories at 
a value of a which is the same in physical units, i.e. for which ayja is the same, does the 
bare coupling vary as 1/iVc) i.e. does (3 = 2Nc/g‘^ oc We perform this comparison for 
+ = 10.9,11.1. (Eor convenience we shall label (3 by the value of Nc, i.e. we write it as 
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(5^^-) To find the corresponding values of (3 in SU{2) and SU{3) we simply interpolate 
between the values provided in (for example) [^]. Doing so we find that the values of f3 
corresponding to (3^ = 10.9,11.1 are Ps ~ 5.972(18), 6.071(24) and j32 2.442(9), 2.485(11) 

respectively. If we simply scale Pi by then what we would have expected to obtain 
is Ps ~ 6.131,6.244 and P 2 — 2.725,2.775 respectively. Superficially the numbers look to 
be in the right ballpark, but in fact the agreement is poor. For example P 2 = 2.725 and 
P 2 = 2.442 correspond to values of a^Ja that differ by about a factor of 3. 

This disagreement should not, however, be taken too seriously, since it is well-known that 
the lattice bare coupling is a very poor perturbative expansion parameter. It is known that 
one can get a much better expansion parameter if one uses instead the mean-field improved 
coupling, g'j, obtained from by dividing it by the average plaquette, < -^TrUp > p^ . 
Defining = 2Nc/g‘j{a) we find that Pi = 10.9,11.1 correspond to Pi = 6.215,6.474 
respectively. Scaling /3| by N‘^ we would expect the equivalent SU (3) and SU (2) couplings 
to be given by P^ = 3.496,3.642 and Pi = 1.554,1.619. What we actually find is that 
the equivalent couplings are Pi ~ 3.527(22), 3.649(28) and Pi = 1.561(10), 1.613(12). The 
agreement is now excellent. That is to say, if the SU{Nc) mean-field improved bare- 
coupling is defined on a length scale that is related to the physical length scale {\/(7) by 
some constant factor, then it varies as g^ oc 1/iVc. This is, of course, the usual diagram- 
based expectation. 

In Fig. § we plot the scalar and tensor glueball masses, in units of as a function of 
Nc- For Nc = 2,3 we have used the continuum values. For W = 4 the calculations are not 
precise enough to permit an extrapolation to the continuum limit and so we simply present 
the values that we obtained at P = 10.9 and 11.1. (We do not use the P = 10.7 values 
since they have large errors and there is the danger that the scalar mass may be reduced 
by its proximity to the critical point at the end of the bulk transition line.) Although the 
W = 4 errors are quite large, it certainly seems that there is little variation with W for 
Nc > 2 and any dependence appears to be consistent with being given by a simple 1/Nc 
correction. The fact that these mass ratios appear to have finite non-zero limits, implies 
that the large-A^c theory is confining. 

As mentioned earlier we have also calculated the topological susceptibility. In Fig. |3| we 
plot the dimensionless ratio a function of Nc- Once again the Nc=2 and 3 

values are continuum extrapolations of lattice values [^, [^, while in the case of 5t/(4) 
we simply display the lattice values obtained at /3=10.9 and 11.1. We remark that for 
SU{A) one expects, semiclassically, very few small instantons and this is confirmed in 
our cooling calculations. This has the advantage that the lattice ambiguities that arise 
when instantons are not much larger than a are reduced as compared to SU{3), and 
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dramatically reduced as compared to SU{2). This implies that the interesting large-A^c 
physics of topology (and the related meson physics) should be straightforward to study. 

4 Conclusions. 

We have calculated the mass spectra and string tensions of gauge theories with W = 2,3,4 
in 3 dimensions. We have found that there is only a small variation with W and this can 
be accurately described by a modest 0{l/N^) correction. That is to say, such theories 
are close to their Nc = oo limit for all values of Nc> 2. We find that the large-W theory 
is confining and that g'^ oc 1/W when expressed in physical units. This confirms, in a 
fully non-perturbative way, expectations arrived at from analyses of Feynman diagrams. 
It simultaneously provides a unified understanding of all our SU (W) theories in terms of 
just the one theory, SU{oo), with modest corrections to it. In practical terms this means 
that, from the parameters in our Tables, we know the corresponding masses for all values 

of W. 

Our calculations in 4 dimensions, while quite preliminary, suggest that the situation is the 
same there as in 3 dimensions. 
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Nc 


2 

0.3350 (15) 

3 

0.5530 (20) 

4 

0.7564 (45) 


Table 1: The D = 2 + 1 SU{Nc) confining string tension. 


/3 

^OO 

Rslope 

CL(%) 

0 ++ 

0.805 ( 13) 

-0.06 ( 8) 

90 

0 ++* 

1.245 ( 27) 

-0.41 (14) 

85 

0 -+ 

1.788 ( 88) 

-0.48 (56) 

40 

2 ++ 

1.333 ( 29) 

-0.08 (18) 

45 

2 “+ 

1.340 ( 40) 

-0.01 (24) 

12 

1 ++ 

1.946 ( 75) 

-0.59 (47) 

95 

1 -+ 

1.919 (115) 

-0.18 (75) 

30 


Table 2; States with C = + va. D = 2 + 1 : R^o 
of the 1/N^ correction in eqn(4). 


lim and Rsiove is the coefficient 

Nc^oo 


G 

^OO 

Rslope 

Q— 

1.18 ( 6) 

0.1 (0.6) 

0 —* 

1.47 (10) 

0.3 (1.1) 

0 +- 

1.98 (28) 

-0.4 (2.7) 

2— 

1.52 (14) 

0.9 (1.4) 

2 +- 

1.58 (13) 

-0.4 (1.3) 

1 — 

1.85 (15) 

-0.3 (1.5) 

1 +- 

1.78 (23) 

1.3 (2.3) 


Table 3; As in Table 2 but for states with C = —. 


10 






/3 

a^a 

amo++ 

am2++ 

10.7 

0.296 (14) 

0.98 (17) 

1.78 (34) 

10.9 

0.229 ( 7) 

0.77 ( 8) 

1.20 (10) 

11.1 

0.196 ( 7) 

0.78 ( 6) 

1.08 (10) 


Table 4: S'C/(4) in 4 dimensions; masses calculated at the values of (5 shown. 
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Figure 1; Some continuum glueball masses, in D = 3, for 2,3,4 colours: 0'’'^(»), 0'’'^*(x), 
2++(*), 0 '■(o) and linear fits. 
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Figure 2; Lightest scalar (•) and tensor (o) glueball masses in Z) = 4. Continuum values 
for iVc = 2,3 and lattice values {(3 = 10.9 and (3 = 11.1) for = 4. 



Figure 3: The topological susceptibility; continuum values for A^c = 2, 3 and lattice values 
{(3 = 10.9 and [3 = 11.1) for Nc = 4. 
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